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On  the  Determination  of  the  Index 
of  Refraction  of  the  Atmosphere 


I.  INTRODUCTION 

It  is  the  intent  of  this  article  to  relate  the  index  of  refraction  of  optical  and 
radio  waves  to  the  distribution  of  atmospheric  species  and  to  demonstrate  that 
with  today's  technology  (using  rocket  measurements  of  density,  temperature,  and 
winds)  it  is  conceivable  that  the  index  of  refraction  in  the  mesosphere  and  strato¬ 
sphere  may  be  directly  measured  under  turbulent  conditions  and  compared  to 

those  values  determined  using  incoherent  radars  and  optical  devices. 

1  2 

As  in  the  Bean  and  Dutton  discussion,  Debye  showed  that  molecular  polari¬ 
zability  is  composed  of  two  effects:  one  due  to  the  distortion  of  all  the  molecules 
by  the  impressed  field  and  the  other  arising  from  an  orientation  effect  exerted 
upon  polar  molecules.  The  polarization  P  of  a  polar  liquid  under  the  influence  of 
a  high-frequency  radio  field  is  given  by  Eq.  (1). 


(Received  for  publication  23  February  1983) 

1.  Bean,  B.R.,  and  Dutton,  E.J.  ( 1966)  Radio  meteorology,  Nat.  Bur.  Stand. 

Monograph  92. 

2.  Debye,  P.  (1957)  Polar  Molecules,  Dover  Publishing  Co.,  New  York. 
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where 

e  is  the  dielectric  constant, 

M  is  the  molecular  weight, 
p  is  the  density  of  the  liquid, 

N  is  Avogadro's  number, 

a  is  the  average  polarizability  of  the  molecules  in  the  liquid  assuming  no 
interaction  between  molecules, 

p  is  the  permanent  dipole  moment, 

k  is  Boltzmann's  constant, 

T  is  the  absolute  temperature, 

t  is  the  relaxation  time  required  for  external  field-induced  orientations  of 
the  molecules  to  return  to  random  distribution  after  the  field  is  removed, 
and 

u  =  2af  where  f  is  the  frequency  of  the  external  field. 

For  frequencies  less  than  100  GHz,  uT  «  1  so  that  Rq.  (1)  may  be  written 


e  -  1  JM  4*N  p2 

e  +  2  p  3  1%  3kTj  • 


(2) 


For  non-polar  gases  (p  =  0),  this  equation  reduces  approximately  to 
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(3) 


Now  from  the  definition  of  M  and  N,  Kq.  (3)  becomes 


e  -  1  s  4*0  n 
o  a 


where  ng  is  the  atmospheric  number  density. 


(4) 
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With  the  inclusion  of  polar  gases  Eq.  (2)  becomes 


(e 


1).  ~  4 *n 
l  a. 


ao  3kT 
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(5) 


where  n  is  the  number  density  of  species  i. 

l 


For  atmospheric  gases,  the  value  of  the  permeability  (fj)  may  be  assumed  to 
be  unity,  with  the  exception  of  water  vapor.  Therefore, 


■ 
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The  equation  for  the  index  or  refraction  n.  s  J  /je  then  equals 


2<n.  -  1)  = 
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or 


N  =  (n.  -  1)106  =  (  ~  1  jx  10G 

2 

"aB2° 

nH20  3kT  ’  (8) 

where  the  summation  now  includes  water  vapor,  and  N  is  now  the  Dielectric  con¬ 
stant. 

To  compare  this  relation  to  measurement,  we  examine  the  dielectric  constant 
at  one  atmosphere  using  the  density  at  ground  level  and  the  dielectric  constant  of 
02  and  Nj.3  Thus, 

(e  -  1)  •  106  =  4*  X  106  (2.  55  X  1019)<~  1. 7  X  10'24)  =  545  . 

3.  Banks,  P.  M. ,  and  Kockarts,  G.  (1973)  Aeronomy,  Academic  Press. 
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3  2*  x  10 


na  *i  + 


This  compares  well  to  the  value  of  575  quoted  in  Table  1  of  Bean  and  Dutton.  * 
Thus,  the  dielectric  constant  N,  as  determined  from  Eq.  (8),  would  be  ~273, 
without  the  contribution  of  water  vapor.  Where  compared  to  the  measured  value 
of  288,  this  result  demonstrates  that  the  index  of  refraction  n.  or  the  dielectric 
constant  N,  as  derived  from  basic  considerations,  are  for  the  most  part  indepen¬ 
dent  of  temperature  and  a  function  of  atmospheric  species,  with  the  exception  of 
water  vapor. 

In  the  lower  troposphere,  water  vapor  is  important  in  the  determination  of  the 
dielectric  constant  through  the  second  term  in  Eq.  (8).  However,  in  the  upper 
troposphere,  stratosphere,  and  mesosphere  water  vapor  is  not  important  to  radio 
wave  or  optical  scattering  by  neutrals  and  Eq.  (8)  becomes 

N  =  (n.  -  1)  •  106  =  (2?r  x  106)  ^  <n  a.)  .  (9) 


2.  THE  EFFECT  OF  TURBULENCE  UPON  OPTICAL  AND  RADIO  WAVE 
SCATTERING 


Given  (n  -  1)  =  C^a  ng  the  fluctuation  component  can  be  defined  as 


where 

rf  is  the  ensemble  average, 
a  s 

The  turbulent  adiabatic  density  fluctuation  is  thus  defined  as 


where 


n 


a 


is  the  potential  density  and 


7  is  the  ratio  of  the  specific  heats. 
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Thus,  it  can  be  shown  that 
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where  Hp  is  the  atmospheric  pressure  scale  height  (=  kT/mg).  Thus,  it  is  easy 
to  show  that  the  Brunt-Vaissalla  frequency  can  be  expressed  as: 


2  JL  =  g.  3© 

UB  ~  9z  ~  0  9z 
n 

a 


(11) 


where  0  is  potential  temperature. 

Given  a  turbulent  region,  the  turbulent  potential  density  " inertial"  spectrum 
may  be  defined  as 


<("a)2>  =  /  ^n  (k)  dk  =  ^1  €‘1/3  cn  /  k'5/3  dk  (12) 

0  3  a  0 

where 

(  )  describe  the  ensemble  average, 

&n  ^  *s  t*le  sca^ar  spectral  energy  density, 
a 

e  is  the  rate  of  viscous  dissipation, 

en  is  the  rate  of  decay  of  the  density  fluctuations, 

k  is  the  wavenumber,  and 

|3j  is  a  constant. 

Similarly,  we  can  define  the  index  of  refraction  fluctuation  spectrum 


<(nl)2>=/  gn.(k)dk  =  P2c'1/3 

0  1 


(13) 
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Making  use  of  Eq.  (8)  we  have  (expressing  all  densities  as  potential  densities 
without  notation) 

(n!2)  =  (Cx  o)2<n^2>  ,  (14) 

where 

contains  all  the  constants. 

We  may  thus  measure  the  amplitude  of  the  turbulent  index  of  refraction  by  the 

direct  measurement  of  atmospheric  density  fluctuations.  This  is  true  for  the 

optical  and  radio  index  of  refraction  in  the  lower  atmosphere  above  the  water  vapor 

region  and  below  ~  55  km  for  the  radio  index  of  refraction.  The  desired  parameter 

to  be  measured  is  the  quantity  C2  ,  the  coefficient  of  the  spectrum,  (=  ft0e  e-1/2), 

n.  ^  n^ 

and  from  Eq.  (14)  we  have  C2  =  (C,  a)2  C2  . 

n.  1  n 

i  a  _ 

Consider  now  the  equation  for  the  production  of  turbulent  heat  flux  (w'61)  where 

w'  is  vertical  turbulent  velocity  fluctuation,  and  O'  is  turbulent  potential  tempera¬ 
ture  fluctuation. 

Ignoring  dissipation,  we  have 


(16) 


where  0  is  the  mean  atmospheric  potential  temperature.  If  the  turbulent  transfer 
and  pressure  terms  are  small  and  in  opposition,  we  have 


<0'2>  S  IHK2) 


In  spectral  form,  we  then  have 


<»'2>  ■  e3<9c‘l/3/  k-5/3  dk  ,  |  ||  „L  e2/3  j  k-5/3  dk 


or  equating  coefficients. 


a  €  *-1/3*0  30  2/3 
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where 


is  a  constant  of  proportionality, 

Cg  is  the  dissipation  rate  of  temperature  fluctuations  by  molecular  diffusion 
2 

((01)  }  is  the  mean  square  temperature  fluctuation,  and 

2  2 
Cq  is  the  representation  of  the  coefficients,  similar  to  C^,  and 

c  .  is  the  constant  of  proportionality  for  the  vertical  component  of  the 

turbulent  velocity. 

2 

Similar  logic  can  be  applied  to  /(n1)  ),  the  neutral  density  fluctuation  giving 


9  n  3n  „ 

W  >*  -f  zT-  <<w,)  > 


and,  for  an  inertial  velocity  spectrum. 


k’5/3dk  • 


Let  us  now  consider  the  component  C  .  This  quantity,  by  definition,  the 
refractivity  structure  constant,  is  given  by'  Tatarskii"*  as  (Balsley  and  Gage'’) 

<  =  5.45  <<n.)2>  l'2/3  ,  (211 

where  Lq  is  the  outer  length  scale  of  the  turbulent  spectra.  Then  the  reflectivity 
(n)  for  backscatter  from  turbulent  fluctuations  by  the  neutral  density  is  (Balsley 
and  Gage'’) 


p  -  0.38  C  X 
ni 


2  > -1/3 


4.  Tatarskii,  V.  I.  (1971)  The  effects  of  the  turbulent  atmosphere  on  wave 

propagation,  U.  S.  Dept.  Commerce,  pp.  74-70. 

5.  Balsley,  B.  B.  ,  and  Gage,  K.  S.  (1980)  The  JUST  ->dar  technique,  Hey.  Pure 

Appl.  Geophys.  118:452. 


where  A  is  the  radar  wavelength.  This  ((Cn  )  )  however  is  nothing  more  than  the 
coefficient  of  the  spectrum  of  the  index  of  refraction  fluctuations,  and  for  a  passive 
neutral  atmosphere,  it  can  be  expressed  as 


and  using  (14)  and  (20) 


n  3n  2/3 

a  a  c  • 


cn.  =  <C1«> 


Substituting  into  Eq.  (22)  we  then  have 


(lf2|  2  na  9na  2/3 
rj  =  (0,  38)  (4ir  )  a  c 


which  directly  relates  the  reflectivity  to  the  dissipation.  Simplifying  the  result 
by  substituting  the  value  of  the  molecular  polarizability  and  constants  (in  the 
Gaussian  mks  system),  we  have 

n.-  =  1.068  X  10‘29  n 
i  a 


Thus, 


c2  =  (1.068  X  10"29)  c2 
"i 


c2  =  (1.068  X  10"29)  aL~  €2/3 


—  -30  3 

where  we  have  used  a  =  1.  7  X  10  m  . 

Tim  question  now  is  how  to  determine  the  neutral  index  of  refraction.  Ob¬ 
viously  one  method  is  by  direct  measurement  of  atmospheric  density  to  fine  enough 
resolution  so  that  we  can  determine  the  local  regions  where  turbulence  is  creating 
these  fine  scale  fluctuations  of  temperature,  velocity,  and  density.  Also  obviously, 
it  would  be  a  prohibitive  venture  in  terms  of  cost  and  time  to  acquire  these  data  to 
the  necessary  accuracy  for  this  project.  Therefore,  to  determine  probable  esti¬ 
mates  of  the  turbulent  contribution  to  the  local  index  of  refraction,  we  must  utilize 
existing  bodies  of  data  that  contain  winds,  temperature,  and  density  measured  to 
a  fairly  fin^  altitude  resolution.  The  only  existing  sequential  data  over  a  large 


enough  area  and  altitude  span  are  the  "Meteorological  Rocket  Network"  and  the 
"Rawinsonde"  data  bases.  However,  given  these  data,  the  next  question  is  how  to 
determine  the  probable  turbulent  contribution  to  the  index  of  refraction,  since 
these  data  in  the  troposphere,  stratosphere  and  mesosphere  have  a  coarse  altitude 
resolution  compared  to  the  turbulent  fluctuation  scale.  There  have  been  two 
approaches  to  this  problem.  These  have  been  formulated  (1)  to  determine  the 

turbulence  contribution  to  the  index  of  refraction  observed  by  radars  (see  the 

5 

reference  list  in  Balsley  and  Gage  )  in  the  troposphere  and  stratosphere,  and 
(2)  to  determine  directly  the  turbulent  parameters  at  local  atmospheric  sites, 
given  the  condition  that  the  local  atmosphere  is  in  a  state  of  destabilization  as 
determined  by  the  calculation  of  the  local  Richardson  number.  It  is  to  this  second 
aspect  that  we  will  now  address  our  efforts. 

r  7 

Zimmerman  and  Murphy,  °  assuming  that  Deacon's  analysis  of  the  boundary 
layer  measurements  of  the  ratio  (w'/U)  (where  w'  is  the  vertical  turbulent  velocity 
and  U  is  the  mean  wind  velocity)  as  a  function  only  of  Richardson  number  could  be 
extrapolated  from  the  boundary  layer  to  the  upper  atmosphere,  determined  the  in- 
situ  turbulent  intensities,  dissipation  rates,  and  diffusivities  that  seemed  quite 
reasonable  when  compared  to  other  measurements  of  these  parameters.  Subse- 

g 

quently,  Zimmerman  and  Keneshea  reasoned  that  the  energy  balance  equation 
would  be  a  more  accurate  representation  of  turbulence  and  probably  would  be  a 
more  altitude  invariant  formulation  than  the  vertical  turbulent  to  horizontal  mean 

7 

velocity  ratio  of  Deacon.  Using  this  reasoning  they  determined  not  only  the  dif¬ 
fusivities  and  dissipation  rates,  but  also  the  local  heat  and  momentum  transfer 

9 

rates.  Comparing  these  analyses  (Figure  1)  to  that  of  Ebel,  who  had  deduced  the 
required  momentum  and  heat  sourees  and  sinks  from  the  Groves*®  model  atmos¬ 
phere,  they  found  quite  good  agreement  of  both  the  phase  and  amplitude  required 
by  Ebel.  In  a  more  recent  work  Zimmerman  and  Keneshea,  **  demonstrate 
* 

This  does  not  imply  that  these  are  the  only  approaches,  since  there  is  a  signif¬ 
icant  body  of  reports  on  the  diffraction  of  optical  wavelength  by  the  atmosphere. 

6.  Zimmerman,  S.  P. ,  and  Murphy,  E.  A.  (1971)  Stratospheric  and  Mesospheric 

Turbulence,  Dynamie  and  Chemical  Coupling,  Reidel  Publishing,  Dordrecht, 
Holland. 

7.  Deacon,  E.  L.  (1959)  Inter.  J,  An  Pollution  2:92. 

8.  Zimmerman,  S.  P. ,  and  Keneshea,  T.J.  ( 198 1)  Turbulent  heating  and  transfer 

in  the  stratosphere  and  upper  mesosphere.  Handbook  for  MAP,  vol  2, 
University  of  Illinois  Press. 

9.  Ebel,  A.  (1974)  Tellus  XXVI.-325. 

10.  Groves,  G.  V.  (1969)  Wind  models  from  60  to  130  km  altitude  for  different 

months  and  latitudes,  J.  Brit.  Interplanet.  Soe.  22:285-307. 

11.  Zimmerman,  S.  P. ,  and  Keneshea,  T.J.  (1983)  Turbulent  heating  and  transfer 

in  the  stratosphere  and  mesosphere  (submitted  to  J.  Atmos.  Terr.  Phys.) 
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FLUX  RICHARDSON  NUMBER 


Figure  1.  The  Normalized,  Non-dimensional  Heat  Flux 
(w'01)  Comparing  the  Data  From  5.  66  m  to  20  km 


that  indeed  a  non-dimensional  heat  flux  and  Reynolds  Stress  may  be  determined 
from  the  boundary  layer  data.  This  non-dimensional  heat  flux  is  shown  to  be  fairly 
invariant  with  altitude  (Figure  2),  at  least  to  20  km.  Thus,  following  that  work  we 
can  determine  the  fluctuation  components  necessary  to  explain  the  local  index  of 
refraction  effects  due  to  turbulence  in  the  stratosphere  and  lower  mesosphere,  if 
we  accept  the  logic  that  we  can  extrapolate  this  non-dimensional  heat  flux  and 

9 

Reynolds  Stress  to  the  higher  altitudes.  The  determinations  of  Cn  using  Eqs.  (21) 
and  (23)  are  that  from  Tatarskii's  analysis  and  that  suggested  here1,  respectively. 

From  Tatarskii 


Cn.  =  5.  45((n!)2)  L~2/3  , 


where  Lq  is  determined  from  the  spectral  equations,  that  is. 


Figure  2.  The  Annual  Average  of  the  Divergence  of  the 
Turbulent  Heat  Flux  Added  to  the  Rate  of  Dissipation 


where 


where  q'  is  turbulent  velocity  and  w',  v'  and  u'  are  the  components.  Thus, 


L  = 


o 


and  ((njF)  is  determined  from 


((n.)2)  --  (Ca)2((n<a 


)2> 


=  (Co) 


2  «e'r>  ~2 

02  na 


(27) 


12.  Groves.  G.  V.  (1971)  Atmospheric  Structure  and  Its  Variations  in  the  Region 
From  25  to  120  km.  AFCRL-G4-0410,  AD  737794. 
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1 


Therefore, 


2  -29  2  (W)  )  ^,2 

C  =  (1.068  X  10  *  )  •  5.45  - -r-  n  • 

ni  a 


Relation  (23)  simply  states 


C2  =  (1.068  X  10'29)  a  -  c2/3 

n.  1  g  8z 


where  e  is  determined  from  the  energy  balance.  Table  1  displays  the  input  condi¬ 
tions  and  results  of  using  the  above  two  relations  on  data  from  the  Meteorologic 
Rocket  Network. 


4 

Table  1.  The  Parameters  for  This  Work  as  Compared  to  Tatarskii's 


This  Work 

Tatarskii 

n  <m'3) 

3n/3z  (m  4) 

<K) 

\soc  / 

C2  (m-2> 

l 

n  (m'3> 

L  (m) 
o 

C2  (nr2) 

l 

1.  2.2324 

6.  1918 

1.  94 

1.  78'18 

721 

1.  922 

l.Sif1" 

ro 

CO 

1.719 

2.  <i'2 

2.  56'17 

1.  962 1 

1.352 

1.7"17 

Note:  The  power  relation  is  10  . 


As  observed,  both  methods  compare  quite  well,  but  Kq.  (23)  is  the  simpler  rela¬ 
tion  to  work  with  and  easier  to  determine. 

Future  work  will  be  to  analyze  the  existing  data  base  to  determine  probability 
2 

factors  of  C  as  functions  of  altitude  and  time  by  turbulence  per  se.  Later  work 

would  be  to  supplement  these  data  with  results  from  analyses  using  the 

13 

Van  Zandt  et  al  formulation  '  for  the  probable  distribution  of  turbulence  given 
measurements  of  Richardson  number  that  indicate  a  non-turbulent  volume. 


13.  Van  Zandt,  T.  E. ,  Gage,  K.  S. ,  and  Warnock  ,  J.  M.  (1981)  An  improved 

model  for  the  calculations  of  profiles  of  C2  and  e  in  the  field  atmosphere 
from  background  profiles  of  wind,  temperature,  and  humidity,  20th 
Conference  on  Radar  Meteorology,  Boston,  Mass. 


References 


1.  Bean,  B.  R.  ,  and  Dutton,  E.J.  (1966)  Radio  meteorology,  Nat.  Bur.  Stand. 

Monograph  92. 

2.  Debye,  P.  (1957)  Polar  Molecules,  Dover  Publishing  Co. ,  New  York. 

3.  Banks,  P.  M. ,  and  Kockarts,  G.  (1973)  Aeronomy,  Academic  Press. 

4.  Tatarskii,  V.  I.  (1971)  The  effects  of  the  turbulent  atmosphere  on  wave 

propagation,  U.  S,  Dept.  Commerce,  pp.  74-76. 

5.  Balsley,  B.B.,  and  Gage,  K.  S.  (1980)  The  MST  radar  technique,  Rev.  Pure 

Appl.  Geophys.  118:452. 

6.  Zimmerman,  S.  P. ,  and  Murphy,  12.  A.  (1971)  Stratospheric  and  Mesospheric 

Turbulence.  Dynamic  and  Chemical  Coupling,  Reidel  Publishing,  Dordrecht, 
Holland. 

7.  Deacon,  E.  L.  (1959)  Inter.  J.  An  Pollution  2:92. 

8.  Zimmerman,  S.  P. ,  and  Keneshea,  T.J.  (1981)  Turbulent  heating  and  transfer 

in  the  stratosphere  and  upper  mesosphere.  Handbook  for  MAP,  vol  2, 
University  of  Illinois  Press. 

9.  Ebel,  A.  (1974)  Tellus  XXV  1:3 2 5. 

10.  Groves,  G.V.  (1969)  Wind  models  from  60  to  130  km  altitude  for  different 

months  and  latitudes,  J.  Brit.  Interplanet.  Soc,  22:285-307. 

11.  Zimmerman,  S.  P. ,  and  Keneshea,  T.J.  (1983)  Turbulent  heating  and  transfer 

in  the  stratosphere  and  mesosphere  (submitted  to  J.  Atmos.  Terr,  Phys.) 

12.  Groves,  G.V.  (1971)  Atmospheric  Structure  and  Its  Variations  in  the  Region 

From  25  to  120  km.  AFCRL-64-0410.  AD  737794. 

13.  Van  Zandt,  T.  E. ,  Gage,  K.S.,  andWarnock,  J.M.  (1981)  An  improved  model 

for  the  calculations  of  profiles  of  C^  and  c  in  the  field  atmosphere  from 
background  profiles  of  wind,  temperature,  and  humidity,  20t.h  Confer  en<  e 
on  Radar  Meteorology,  Boston,  Mass. 


17 


